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ABSTRACT 

The validity of the stratospheric aerosol measurements made by the satellite sensors SAM II 
and SAGE has been tested by comparing their results with each other and with results 
obtained by other techniques (lidar, dustsonde, filter, impactor) . The latter type of com- 
parison has required the development of special techniques that (1) convert the quantity 
measured by the correlative sensor (e.g. particle backscatter, number, or mass) to that 
measured by the satellite sensor (extinction), and (2) quantitatively estimate the uncer- 
tainty in the conversion process. The results of both types of comparisons show agreement 
within the measurement and conversion uncertainties. Moreover, the satellite uncertainty 
is small compared to aerosol natural variability (caused by seasonal changes, volcanoes, 
sudden warmings, vortex structure, etc.). Hence, we conclude that the satellite measurements 
are valid. 


INTRODUCTION 

The satellite sensors Stratospheric Aerosol Measurement (SAM II) and Stratospheric Aerosol 
and Gas Experiment (SAGE) have been developing a global data base on stratospheric aerosol 
behavior since they were launched in October 1978 and February 1979, respectively / 1—6/ . To 
test the validity of this data base, an extensive program of comparisons has been conducted, 
both between the two satellite sensors and between other sensors and the satellite sensors. 

COMPARISONS BETWEEN SAM II AND SAGE 

As described in more detail by Yue et al. / 7 / , the SAM II measurements are confined to the 
two latitude bands: 64“-80°N and 64°-80°S, whereas the SAGE measurements range from 72°S- 

to 72”N. Thus, there are occasions when the measurements nearly overlap, permitting 
comparisons . 

Although SAM II and SAGE use the same measurement principle (sun photometry through the 
Earth's atmospheric limb, yielding vertical profiles of aerosol extinction), there are sig- 
nificant differences between the instruments used and data-reductlon procedures. For 
example, the sensors fly on separate satellites in different orbits, leading to considerably 
different viewing geometries and sampling durations even when the sampled air volumes are 
adjacent. Also, the SAGE instrument uses a holographic grating to disperse the sunlight, 
while the SAM II instrument uses an interference filter for spectral discrimination. SAGE 
measurements are quantized to 12 bits accuracy (1 part in 4096) , while SAM II measurements 
are quantized to 10 bits accuracy (1 part in 102**) . The higher quantization error for 
SAM II measurements is apparent when the measured aerosol extinction levels are low. 

The handling and processing of the data from each instrument are also different. The 
telemetry and raw data processing for the two satellite systems are separately operated 
under different organizations at NASA's Goddard Space Flight Center. In addition, the data 
reduction schemes for the two experiments are considerably different. For example, the 
geometric height determination for the SAM II-inferred vertical profiles is based primarily 
on the refracted sunshape fitting scheme /8/. For the SAGE data, it is determined by match- 
ing the measured slant-path optical depth profiles at the short wavelength channel to a 
computed profile from the given vertical temperature and pressure profiles. The inversion 
algorithms for the two instruments are also different. SAGE inversion uses a nonlinear 
iterative inversion algorithm where all the four wavelength data are cyclically updated, 
while the SAM II inversion algorithm is basically a linear inverse matrix-type inversion 
applied to the single wavelength measurement. The comparisons between SAM II and SAGE 
1.0-um wavelength inverted data, therefore, can provide a critical assessment of 

05 ' 



2 


the end-to-end performance of Che two measurement systems, and can validate the consistency 
of the inferred aerosol vertical extinction profiles. 

By the end of 1979 there were five occasions when the latitude coverages of SAM II and SAGE 
overlapped. We will show representative comparisons, both between individual profiles and 
between height-vs. -longitude cross sections generated by combining many profiles. 

SAM II and SAGE Individual Profile Comparisons 


Figure 1(a) shows a typical SAM II/SACE comparison for a case of background aerosol concen- 
trations, where the measurement coordinates differed by 0.26° latitude, 0.86” longitude, and 
0.06 hours. (The coordinates give the location and time when each satellite's Sun-viewing 
limb path was tangent at 20 km altitude.) The error bars shown are ilo, derxved from the 
error analysis of Chu and McCormick /8/. They include uncertainties in pointing, radio- 
metric imprecision, digitizer truncation, and uncertainties m the molecular density profile 
at the measurement location. The profiles are typical in that they show a steep decrease of 
aerosol extinction near the tropopause height (indicated by the arrow), followed by a region 
of nearly constant extinction extending from about 2 km above the tropopause to about 8 km, 
followed by a relatively steep drop-off. 

As can be seen, within the stratosphere the SAM II and SAGE profiles agree to within error 
bars. Below the tropopause, differences often exceed error bars. We attribute this behavior 
(which is typical of other profile comparisons) to the patchiness of the tropospheric aerosol 
structure, coupled with the fact that SAM II and SAGE tropospheric limb paths can be sepa- 
rated by considerably more chan the differences stated above and shown in Figure 1. (Recall 
chat those differences apply to the 20-km tangent height.) In cases where tropospheric 
clouds are present, SAM II/SAGE differences can greatly exceed those shown in Figure 1(a) 
because of the large inhomogeneities of clouds. For other examples and further discussion, 
see / 7/ . 

Figure 1(b) shows a SAM II/SAGE comparison for a case of unusual background stracospheric 
aerosol structure. Note that, in this case, particulate extinction decreases rapidly from 
about 2 * 10 -1 * km” 1 at 11 km altitude to about 2 * 10“ s at 18 km altitude, rather than 
staying nearly constant which is far more typical behavior in this height region (cf. Fig- 
ure 1(a)). This interesting phenomenon was not a single example but continued for several 
days. The mechanism behind this phenomenon is still under investigation /6/ ; however, it 
is possible chat this is the result of a strong polar vortex and is typical of the air mass 
inside che vortex experiencing a subsidence and sedimentation. The Important point here is 
that SAM II and SAGE are in excellent agreement in documenting this atmospheric behavior. 



AEROSOL EXTINCTION l/Vm 



AEROSOL EXTINCTION 1/km 


Fig. 1. Selected comparisons between SAM II and SAGE aerosol extinction profile 

measurements . 
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SAM II and SAGE Concour Plot Comparisons 

Figure 2 shows two examples from SAM II and SAGE in which height vs. longitude cross sections 
are compared. For each sensor, these cross sections are generated by inputting to a contour 
plotting routine all the vertical profiles measured in a single day. Numbers shown on the 
contours are in units of 10 -5 km -1 and are separated by multiples of 10° ‘ 2 (about 1.58). As 
in the individual profiles shown in Figure 1, only values of extinction greater than 
10” 5 km~ l are shown; the solid lines represent data obtained by SAGE, and the dashed lines 
are for SAM II. The longitudinal positions where SAGE and Sam II data were taken are indi- 
cated by tick marks at the top and bottom of each figure, respectively. The latitudinal cov- 
erage for each sensor for that day is also listed. The positions of each Measurement indi- 
cated in Figures 2(a) and (b) are those corresponding to a limb path at 20 km altitude. On 
a given day, 15 is about the maximum number of sunrise or sunset measurement opportunities 
for each sensor. 


SAGE-172 63 N-73 14 N) 

AUG 6 1979 

164 5 -88 9 -38 5 119 94 9 145 2 



SAM 11-172 35 N-72 55 Nl 
l»> LONGITUDE deq 
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Fig. 2. Selected comparisons between height/ longitude cross sections of extinction derived 
from SAM II and SAGE measurements. Solid lines: SAGE; dashed: SAM II. 
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The reader should be cautioned that computer-generated contours such as those in Figure 2 
can lead to apparent differences and possible misinterpretations, because profiles at dif- 
ferent longitudes are used in their independent construction. Many times and for various 
reasons, a profile is missing; thus, instead of a 24“ (SAGE) or 26“ (SAM XI) separation 
between profiles for interpolation, twice this longitude interval is used. Sometimes two 
consecutive profiles are missing, causing an even larger interpolation distance. If, for 
example, there is a localized enhancement in the SAM II data that coincides with a data-void 
area for SAGE, a very different contour would be constructed for the two sensors. 

From the plotted contours in Figure 2(a) which represents a very homogeneous condition, it 
can be seen that there is very good agreement between the SAto II and SAGE measurements for 
the altitude region more than 2 km above the tropopause height (indicated by the circles for 
SAGE data and squares for SAM II data) . Altitude differences are in general less chan 1 km 
between contours of the same aerosol extinction, but obtained from different sensors. As 
mentioned earlier, the small differences between SAM II and SAGE measurements of aerosol 
extinction in the region around the tropopause may be due to variable tropopause heights, 
and the variable and patchy clouds associated with them. The small differences in the 
region above 24 km are caused by relatively large systematic errors associated with data in 
that region. Even though there are small differences between SAM II and SAGE results in 
those two regions, they are still in good agreement within their respective uncertainties. 
Most of the differences seen in Figure 2(a) are due to data voids and their effects on the 
contour plotting. 

The aerosol extinction contour plots shown in Figure 2(b) represent the case in which there 
are two regions where aerosol extinction is very different. As described previously, this 
phenomenon is associated with a well-established polar vortex / 6/ . Figure 2(b) vividly 
shows the distortion generated by having data gaps at a particular longitude. Even here, 
however, both contours show the major features of the stratospheric aerosol on that day. 

COMPARISONS OF OTHER SENSORS WITH SAM II AND SAGE 

Of all the nonsatellite aerosol sensors and samplers available to validate SAM II and SAGE 
(e.g. Ildars, optical particle counters, filter collectors, and various types of impactors) , 
none measures extinction, which is the quantity measured by SAM II and SAGE. Instead, the 
nonsatellite sensors measure such quantities as particle backscatter, number, and mass. 

Thus, before using these correlative sensors to validate the satellite sensors, procedures 
must be developed that (1) convert the correlative-measured property to the satellite- 
measured one, and (2) quantitatively estimate the uncertainty in the conversion process. 

Conversion Procedures and Uncertainties 


It i3 well known that there are no unique factors (or conversion ratios) relating extinction 
to mass, number, or backscatter for the stratospheric aerosol because the factors (ratios) 
depend on particle size distribution, composition, and shape — all properties which are known 
to vary. Nevertheless, one can determine the range of these properties by using concurrent 
measurements, physical principles, and climatological information. Having defined the 
ranges of size distribution, composition, and shape, one can then define the corresponding 
ranges of possible conversion ratios. These ranges may be expressed as mean values and 
uncertainties, each as functions of height. 

Figure 3 illustrates schematically our technique for accomplishing and automating the above 
process. The technique is explained In detail in /9-11/. In brief, it consists of the 
following steps: 

1. Define the range of aerosol refraction Indices expected in specified height regions. 

2. Parameterize a wide range of model relative-size distributions in terms of the "channel 
ratio," No.u/Nu. 25 . where N x is the number of particles with radii greater than x ym. 
(We chose channel ratio as the key size-distribution parameter because it is measured by 
dustsondes in most of our correlative experiments, and when measurements are not avail- 
able, there is a large data base from which channel ratio can be estimated. Moreover, 
for the stratospheric aerosol, knowledge of the channel ratio provides a very useful 
constraint on conversion ratio values.) 

3. Perform Mie-scattering calculations for a very large number of combinations of refractive 
indices and size distributions, and organize the results to give mean and standard devia- 
tion curves of conversion ratios (e.g. extinction-Co-mass) as a function of No.is/No.25 
(see /9,10/ for examples). The mean curves give the most probable conversion ratio for a 
particular No.is/No.ts value, and the standard deviations give the la uncertainty 

(i.e. expected rms error) that results from imperfect knowledge of refractive index and 
size-distribution details (see the example in Figure 3(a)). 
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Fig. 3. Schematic of the procedure used to convert measurements of correlative properties 
(e.g., mass, number, backscatter) to estimates of extinction. 

4 . Combine curves of conversion ratio versus No.is/No.25 with measured (or estimated) verti- 
cal profiles of No.is/No .25 (Figure 3(b)) to obtain vertical profiles of conversion ratio 
and its uncertainty (Figure 3(c)). 

5. Multiply measured profiles of the correlative property (number, backscatter, mass; 

e.g. Figure 3(d)) by the corresponding profile of conversion ratio (extinction-to-number , 
extinction-to-mass , or extinction-to-backscatter) , and combine measurement and conversion 
uncertainties, using conventional error-propagation methods (see /9— 10/) . The result is 
the correlative extinction profile exemplified by Figure 3(e). 

Results: Lidar vs. SAM IX; Vol-canlc Aerosol 


Figure 4 shows two examples of comparisons between extinction profiles measured by SAM II 
and those estimated from airborne lidar measurements of backscatter. In both cases, the 
stratospheric aerosol concentration is greatly enhanced over the background levels of 
Figures 1 and 2, primarily because of the El Chichon volcanic injection. Due to this vol- 
canic perturbation, the conversion procedure described above and in Figure 3 was modified 
slightly, both because (1) simultaneous duscsonde measurements of No.is/No.25 were not avail- 
able (and no "typical El Chichon" No 15 /N 0.25 profiles could be estimated from the sparse 
data available elsewhere), and (2) the preponderance of large particles in the El Chichon 
cloud decreases the usefulness of No.is/No .25 in constraining size distribution. Instead of 
relying exclusively on No.is/No.25, 6-channel dustsonde measurements of post-El Chichon size 
distributions were used to define mean and standard-deviation extinction-to-backscatter 
ratios. These were then combined with the lidar-measured backscatter profiles to yield the 
dashed extinction curves and error bars in Figure 4. It can be seen that, in spite of the 
uncertainty in particle size distribution, the resulting uncertainty in lidar- inf erred 
extinction is relatively small (i20Z in the aerosol peak). Within this uncertainty, there 
is a very good agreement between the satellite-measured and lidar-inferred extinction 
profiles. 

Results: Lidar and Dustsonde vs. SAM II, Background Aerosol 

Figure 5 shows results from a validation experiment /10/ conducted at Sondrestrom, Greenland, 
In November 1978, during "background" (nonvolcanic) aerosol conditions. In this case, both 
a lidar and a dustsonde participated, and the conversion procedure shown in Figure 3 was 
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followed exactly. The error bars on the lidar- and dustsonde-inferred extinction curves 
give the combined result of measurement and conversion uncertainty. (For the dustsonde, 
conversion errors dominate the total, for the lidar, measurement errors dominate. See / 10/ 
for details.) For both days (22, 23 Nov.) at all heights above the tropopause, the extinc- 
tion profiles measured by SAM II agree with those inferred from the lidar and dustsonde 
data to within the respective measurement and conversion uncertainties. Below the tropo- 
pause on 23 Nov. 1978 (Figure 5) , the SAM II-measured extinction exceeds chat measured by 
the dustsonde. We believe this is caused by the SAM II limb path intercepting patchy cirrus 
clouds which were not efficiently sampled by the dustsonde. 


Results: Lidar. Dustsonde, Filter, and Impactor vs. SAM II and SAGE. Background Aerosol 


Figure 6 shows the results from a validation experiment /11,12/ conducted at Poker Flat, 
Alaska in July 1979, also during background aerosol conditions. The procedure illustrated 
in Figure 3 was used to derive extinction profiles and uncertainties from the lidar, dust- 
sonde, and filter measurements of backscatter, number, and mass, respectively. We have also 
derived extinction values from the size distributions measured by wire impactor. However, 
our error analysis for these values differs from that for the other sensors because of the 
more detailed nature of the vire-impactor data. Specifically, the wire impactor measures 
differential particle number densities (dN r /dr) at a large number (-15) of radius intervals, 
instead of directly measuring integral properties (such as Ng.u, N0.25, mass, and back- 
scatter). Thus, the wire impactor yields detailed size-distribution information, but each 
dN r /dr data point has associated uncertainties, amounting to i8Z to t20Z in r and highly 
variable amounts in dN r . Because of the steepness of background stratospheric partlcle-sizf 
distributions, the uncertainties in r lead to considerably larger uncertainties in some 
wire-impactor-derived integral properties (such as Ng.is/Ng.2s, mass, and extinction). 


A complete analysis of the la uncertainties in these impactor-derived integral properties is 
very complicated; however, to illustrate the magnitudes possible, we show in Figure 6 two 
wire-impactor-inferred extinction values at each measurement height. One value is derived 
from the "best-estimate" impactor- inferred size distribution, and the other is derived from 
this same distribution with all radii reduced by 15Z. Note chat this small radius shift 
(which is within the wire— impactor measurement uncertainty — 95Z confidence limits), reduces 
extinction by a factor of 2 to 3. Thus, the more detailed wire-impactor data, when used as 
an indirect means to estimate extinction (and certain other integral properties) , yield 
relatively large uncertainties. 


Figure 6 also shows the profiles of extinction measured by SAM II and SAGE. To within the 
associated measurement and conversion uncertainties, the SAM II and SAGE profiles agree with 
each ocher and with the values derived from the simultaneous dustsonde, filter, wire- 
impactor, and lidar measurements. 



EXTINCTION COEFFICIENT - km' 1 


Fig. 6. Comparison of extinction profiles measured by SAM II and SAGE with those derived 
from dustsonde (DS) , filter, lidar, and Ames wire-impactor (AWI) measurements. For each 
horizontal pair of AWI data points, the point on Che left results from reducing all parti- 
cle radii by 15Z. 
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CONCLUSIONS 

The SAM II and SAGE measurements agree with each other and with values inferred by other 
techniques to within the respective measurement and conversion uncertainties. Moreover, the 
satellite measurement uncertainty is small compared Co aerosol natural variability. Hence, 
we conclude that the satellite measurements are valid. 
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